arXiv: 1504.04265vl [cond-mat.mes-hall] 16 Apr 2015 


Thermoelectric and thermospin transport in a ballistic junction of graphene 
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We consider theoretically a wide graphene ribbon, that on both ends is attached to electronic 
reservoirs which generally have different temperatures. The graphene ribbon is assumed to be 
deposited on a substrate, that leads to a spin-orbit coupling of Rashba type. We calculate the 
thermally induced charge current in the ballistic transport regime as well as the thermoelectric 
voltage (Seebeck effect). Apart from this, we also consider thermally induced spin current and spin 
polarization of the graphene ribbon. The spin currents are shown to have generally two components; 
one parallel to the temperature gradient and the other one perpendicular to this gradient. The latter 
corresponds to the spin current due to the spin Nernst effect. Additionally, we also consider the 
heat current between the reservoirs due to transfer of electrons. 

PACS numbers: 72.25.Fe, 78.67.Wj, 81.05.ue, 85.75.-d 


I. INTRODUCTION 

Low energy electronic states in graphene - a two- 
dimensional hexagonal lattice of carbon atoms ~ are usu¬ 
ally described by the relativistic Dirac model>i Unique 
transport properties of graphene, especially the high elec¬ 
tron velocity and tendency to avoid electron scatter¬ 
ing due to the Klein effect, make graphene an excel¬ 
lent material for future applications in nanoelectronicsi^ 
These properties of graphene also facilitate practical re¬ 
alization of ballistic junctions with graphene4i^ Indeed, 
such junctions have been extensively studied in recent 
yearsi^— It has been shown, for instance, that junc¬ 
tions of two graphene parts corresponding to different 
Rashba coupling parameters exhibit interesting transport 
featuresi^ When the graphene is additionally magnetized, 
e.g. due to coupling to an insulating magnetic substrate, 
a large anisotropic magnetoresistance effect can be then 
observed 

Remarkably less theoretical and experimental work has 
been done up to now on thermoelectric properties of 
graphene, though interest in these properties is growing 
recentlyjiic— Theoretical works were focused mainly on 
the diffusion transport regime in graphene with impu¬ 
rities and other structural defects. It has been shown, 
for instance, that at certain conditions the Wiedemann- 
Franz law can be violated in graphenei^ Moreover, reso¬ 
nant scattering from impurities with short-range poten¬ 
tial may lead to an enhanced Seebeck coefficient, when 
the chemical potential is in the neighborhood of the 
resonancesOf particular interest are currently ther¬ 
moelectric properties of graphene nanoribbons, which 
may exhibit an enhanced thermoelectric efficiencyi^^— 
This efficiency may be additionally enhanced by certain 


structural defects like antidots, for instancei^^ In addi¬ 
tion, there is currently a great interest in spin related 
thermoelectric phenomena in nanoscale systems, includ¬ 
ing also graphene nanostructures. It has been shown, 
among others, that graphene nanoribbons with zigzag 
edges can exhibit not only conventional but also spin 
thermoelectricityji^ The latter corresponds to a spin volt¬ 
age generated by a temperature gradient. Various phys¬ 
ical phenomena associated with thermally induced spin 
and heat currents are also of current interest 

One may observe recently an increasing interest in the 
ballistic transport regime j^^i^^ This is due to the possibil¬ 
ity of a long mean free path t in 2D electron gas^^, where 
£ ~ 3 /rm has been already reached,— and in graphene 
nanoribbons, where £ > 10 /rm has been reported— The 
main objective of this paper is a theoretically descrip¬ 
tion of thermoelectric and thermospin transport prop¬ 
erties of ballistic graphene junctions. We calculate the 
thermoelectrically induced charge and spin currents in a 
graphene ribbon of length L, which is attached to two 
electronic reservoirs. The ribbon length is assumed to be 
smaller than the corresponding electron mean free path 
t. The main focus in the paper is on the spin effects 
due to spin-orbit interaction in graphene. Since the in¬ 
trinsic spin-orbit coupling in graphene is very small, it 
is neglected here. In turn, the Rashba spin-orbit cou¬ 
pling related to the influence of a substrate can be rela¬ 
tively strong, and therefore it is included in our consid¬ 
erations. In the presence of Rashba spin-orbit coupling 
in graphene, temperature gradient can generate not only 
the spin current (spin Nernst effect) but also a spin den¬ 
sity. It is well known, that spin current can be then gener¬ 
ated also by an electric field (spin Hall effect). Similarly, 
the spin polarization may be induced by the electric held 
as well— 
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In section 2 we describe the theoretical model. Nu¬ 
merical results on the thermally induced charge and heat 
currents are presented and discussed in section 3. In turn, 
thermally induced spin current and spin polarization of 
graphene is considered in section 4. Summary and final 
conclusions are in section 5. 


II. MODEL 

We assume the relativistic Hamiltonian for electrons 
in graphene with Rashba spin-orbit coupling, 

H^ = hvFT ■'k + a{(Tj;Ty-ayT^), ( 1 ) 

where vf is the electron velocity in graphene, k = 
{kx^ ky) is a two-dimensional wavevector, a is the Rashba 
coupling parameter, while r and cr are the vectors of 
Pauli matrices defined in the sublattice and spin spaces, 
respectively. Hamiltonian (1) describes low-energy elec¬ 
tron states in the vicinity of the Dirac point RT of the 
corresponding Brillouin zone. Hamiltonian for the sec¬ 
ond non-equivalent Dirac point, K', can be obtained from 
Eq.(l) by reversing sign of the wavevector component k^. 


Ti Ti > T2 T2 



X * 


Figure 1: Schematic of a ballistic junction consisting of a 
graphene ribbon of length L and two 2D electronic reservoirs. 
The reservoirs have generally different temperatures, as indi¬ 
cated. 

The electronic band structure described by the Hamil¬ 
tonian (1) consists of four energy bands, 

Ekn = ±a± (ffvpk'^ + ( 2 ) 

where n = 1 — 4 is the band index, with n = 1 (n = 4) 
corresponding to the band of the lowest (highest) energy. 
Each of the bands is parabolic at small wavevectors, k <C 
a/hvF, and has almost linear dispersion for k ^ a/hvF- 
Two of these bands (n = 2, 3) touch at fc = 0, while two 
others (n = 1,4) are separated by an energy gap of width 
equal to 2q:. 

We assume the axis x is along the graphene ribbon of 
length T, as shown schematically in Fig. 1. The ribbon 
is wide enough to neglect size quantization. In turn, the 
length L is smaller than the mean free path i, L ^ £, so 
electronic transport can be considered as fully ballistic. 


For simplicity, we neglect any scattering of electrons in¬ 
side the ribbon. Additionally, we assume the graphene 
ribbon is in contact at the ends with two 2D electronic 
reservoirs, which generally have different temperatures, 
Ti and T 2 , as indicated in Fig.l. Electrons in these 
reservoirs are described by the corresponding equilibrium 
Fermi-Dirac distribution functions, 

fiA^kn) = [exp((ekn - fJ-l,2)/kBTi^2) + 1]“\ (3) 

where fii and fj ,2 are the chemical potentials. Though we 
are focused mainly on thermal effects, we assume that 
and fi 2 can be different in a general situation. Thus, the 
electron system in the ballistic region can be described by 
the distribution functions /(^(ekn) and /^(ckn) for elec¬ 
trons moving from left to the right and from the right to 
left, respectively. In the following we use this distribu¬ 
tion to calculate transport and thermoelectric properties 
of the graphene ribbon, assuming purely ballistic regime. 


III. THERMALLY INDUCED CHARGE AND 
HEAT CURRENTS 

Assume equal chemical potentials in the two electronic 
reservoirs, /ii = /r 2 = /i, and different temperatures, 
Ti > T 2 . The former condition will be relaxed only 
when necessary. Below we calculate charge current due to 
the temperature difference (gradient), and also the cor¬ 
responding electronic contribution to the heat current. 


A. Charge current 

The charge current in the ballistic transport regime, 
flowing along the axis x due to the difference AT = Ti — 
T 2 in the reservoir temperatures, can be calculated with 
the formula 

j = e^^ {kn\vx\kn)[f^{ekn)-f2i£kn)], (4) 

n k 

where e is the electron charge, Vx = vfTx is the elec¬ 
tron velocity operator, and the summation over the 
wavevector k is restricted to the angles, for which the 
x-component of electron velocity, Vxn = (kn|f)x|kn), is 
positive. 

For definiteness, we assume that the chemical potential 
of electrons is positive, ^ > 0. It is clear that due to the 
electron-hole symmetry, the results for /i < 0 can differ 
only in sign from those for ^ < 0. This is because the 
electron velocity Vxn is positive for kx > 0 'u\ the energy 
bands with £kn > 0 and negative for the bands with 
S^kn ^ d- 

The thermoelectric current calculated as a function of 
AT is presented in Fig. 2a for different values of the 
Rashba coupling constant a. Using Eq. (3) one can 
show that the dependence of current on AT is linear for 
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Figure 2: Thermoelectric current in the ballistic regime as 
a function of AT (a) and chemical potential fi (b) for indi¬ 
cated values of the Rashba coupling parameter a, and for 
fj. = 5 meV (a) and AT = 14 K (b). 


a = 0 and AT ^ Ti, /r/Zcs, which is related to the lin¬ 
earity of the density of states in graphene, ^{e) |e|. 

It should be noted that in the case of ordinary 2D elec¬ 
tron gas with parabolic energy spectrum, the thermo¬ 
electric current at these conditions is absent since the 
corresponding 2D density of states, V 2 d, is independent 
of the electron energy, and thus the currents due to par¬ 
ticles and holes compensate each other, so the net cur¬ 
rent vanishes. Thus, the nonzero thermoelectric ballistic 
current in graphene is related to the relativistic energy 
spectrum of this 2D material. As follows from Fig.2a, the 
dependence of the thermoelectric current on the Rashba 
spin-orbit constant a is non-monotonic, and the current 
is maximal when a ~ 

Variation of the thermoelectric current with the chem¬ 
ical potential fj, is shown in Fig.2b. First, the current 
vanishes when the Fermi level is at the particle-hole sym¬ 
metry point, ^ = 0. Second, there is a small peak for 
larger values of a. Position of this peak corresponds to 
the onset of the contribution from another (higher in en¬ 
ergy) band, whose bottom band-edge is at the energy 
s = a. At larger values of the electric current satu¬ 
rates at a constant value (independent of a), since the 
energy spectrum is linear in this region, like in graphene 
without spin-orbit coupling. 

To compare the results for ballistic and diffusive junc¬ 
tions, we note that Eq. (3) also gives the thermoelectric 
current in the diffusive regime if we take AT = U^T and 



Figure 3: Comparison of the thermoelectrically induced cur¬ 
rent in the ballistic and in the diffusive transport regimes, 
calculated for jj, = 5 meV. In the diffusive regime, the electric 
current is proportional to the temperature gradient AT. 



Figure 4: Thermoelectric voltage as a function of AT, cal¬ 
culated for fj, = 5 meV and indicated values of the Rashba 
parameter a. 


substitute L by the mean free path i, as it should be to 
match the ballistic and diffusive results. Thus, if we keep 
VT = const and reduce £, i.e., if we go to the diffusive 
regime by increasing the density of impurities, £< L, we 
decrease the current. The difference between the ther¬ 
moelectric currents for ballistic and diffusive transport 
regimes is shown in Fig. 3. While in the diffusive regime 
the thermo-current increases linearly with the tempera¬ 
ture difference AT, its increase with AT in the ballistic 
limit is faster and nonlinear. 

One can also calculate the thermopower, or equiva¬ 
lently the thermally induced voltage U between the reser¬ 
voirs under the condition of zero charge current, j = 0. 
To determine the voltage U, one can use Eq. (3) with the 
distribution functions corresponding to different electro¬ 
chemical potentials, fii ^ fj, 2 , as in Eq. 2. Then, the 
voltage U can be determined from the condition that the 
thermally induced current is fully compensated by the 
held-induced current. The thermoelectric voltage is then 
given as U = (^i — T 2 )/e. Results of the correspond¬ 
ing numerical calculation are presented in Eig. 4. As 
follows from this figure, the thermoelectric voltage de¬ 
pends on the temperature difference in a nonlinear way, 
whereas the dependence on the Rashba parameter is non- 




































4 



AT (K) 



Figure 5: Heat flux transmitted by electrons for indicated 
values of the Rashba parameter a, calculated as a function of 
AT for /i = 5 meV (a), and as a function of for AT = 14 K 
(b). 

monotonous. 


B. Heat current 

Employing the same method as that applied above in 
the calculation of thermoelectric charge current, one can 
also calculate the heat current associated with transfer of 
electrons between the two reservoirs. The corresponding 
formula for the heat flux from the left reservoir (of tem¬ 
perature Ti) to the right one (of temperature T 2 ) can be 
written in the form 

JQ = 1 EE' (kn|{(i7k - /r), fla;}|kn) 

n k 

X [fi (Ekn) - f2 (^kn)] , (5) 

where {A^ B} = Al3 + BA for any two operators A and 
13. 

Dependence of the heat flux Jq on the temperature 
difference AT is presented in Fig. 5a for different val¬ 
ues of the Rashba parameter a. Similarly as in the 
case of charge current, the heat current increases non- 
linearly with AT, and also depends nonmonotonously on 
the Rashba parameter a, compare Fig.2a and Fig.5a. 

Figure 5b in, turn, shows the dependence of the heat 
current Jq on the chemical potential /i. As follows from 




Figure 6: Thermoelectrically induced spin polarization as a 
function of AT (a) and as a function of p (b), calculated for 
indicated values of the Rashba parameter and for fj, = 5 meV. 
(a) and AT = 14 K (b). 


this figure, the dependence on ^ is linear at large values 
of contrary to the behavior of charge current which 
saturates at large /i (see Fig.2b). This difference appears 
because the contribution to heat current from a trans¬ 
ferred electron depends on its energy, e ^ while the 
corresponding contribution to charge current is indepen¬ 
dent on this energy. 


IV. THERMALLY INDUCED SPIN 
POLARIZATION AND SPIN CURRENT 

As in the preceding section, we assume equal chemical 
potentials in the two electronic reservoirs, /ii = /i 2 = T- 
Below we calculate spin polarization of electrons and spin 
current, both induced by a temperature difference AT. 


A. Spin polarization 

It is already well known that spin-orbit interaction in 
the presence of either electric field or temperature gradi¬ 
ent can induce spin polarization of conduction electrons. 
This effect has been studied theoretically for the usual 2D 
electron gas with parabolic energy spectrum and Rashba 
spin-orbit coupling^^— as well as in graphene^i In the 
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ballistic junction considered in this paper, the spin den¬ 
sity can be calculated using the formula 

IM [fl (ekn) - f2 (^kn)] . (6) 

n k 

The corresponding numerical results are presented in 
Fig. 6, where the thermally-induced spin polarization Sy 
is shown as a function of AT (Fig.6a) and as a function 
of /i (Fig.6b). The induced spin polarization is in the 
graphene plane, and is normal to the temperature gradi¬ 
ent, similarly as in the case of 2D electron gas4^ Physi¬ 
cal mechanism of the spin polarization in graphene with 
Rashba spin-orbit coupling is also similar to the mech¬ 
anism of spin polarization in 2D electron gas. Indeed, 
there is a nonzero spin polarization of an electron in the 
eigenstate |kn), which is perpendicular to the wavevec- 
tor k. The magnitude of this spin polarization is small 
for k <C ajhvF and is equal to its maximum value equal 
to h/2 for k ^ a/hvp^ A nonzero spin polarization 
Sy appears due to the imbalance of the distribution of 
electrons with kx > 0 and fca, < 0.— 


B. Spin currents 

The temperature difference between the electron reser¬ 
voirs can also generate a spin current flowing parallel 



AT (K) 



Figure 7: Thermally induced spin currents (a) and Jjj 
for different values of a, calculated as a function of AT for 
fj. = 5 meV. (a). 



/Li (meV) 


Figure 8: Thermally induced spin currents (a) and Jy 
for different values of a, calculated as a function of for 
AT = 14 K. 


to the temperature gradient as well the spin current Jy 
flowing perpendicularly to the gradient. Here the upper 
index indicated the spin component associated with the 
spin current, while the lower index indicated the orien¬ 
tation of the spin current flow. Both components of the 
spin current can be calculated using the same method as 
in the case of charge current. The relevant formula takes 
now the form 

(kn|{crc,D}|kn) [fliSkn) - /2‘^(£kn)]-(7) 

n k 

The numerical results for both and Jy calculated 
as a function of the temperature difference AT are pre¬ 
sented in Fig. 7 and as a function of the chemical poten¬ 
tial /i in Fig. 8. The mechanism of a nonzero component 
J]^ is related to the spin polarization of electrons due to 
the temperature gradient and Rashba spin-orbit interac¬ 
tion, as calculated and discussed above. These spin po¬ 
larized electrons are transferred between the two electron 
reservoirs, giving rise to the spin current J^. In turn, the 
other spin current component, J®, corresponds to the 
thermally induced spin Hall effect, called also the spin 
Nernst effect. This effect consists in a spin current gen¬ 
eration by a temperature gradient. The induced current 
flows then perpendicularly to the temperature gradient. 
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V. SUMMARY 

We have analyzed thermoelectric and thermospin ef¬ 
fects in a ballistic graphene ribbon attached to two elec¬ 
tronic reservoirs of different temperatures. The graphene 
ribbon was assumed to be deposited an a substrate 
that generated a strong spin-orbit coupling of Rashba 
type. We have calculated not only the thermally induced 
charge current between the two reservoirs, and the asso¬ 
ciated thermoelectric voltage, but also thermally induced 
spin polarization and spin current. Numerical results on 
the charge current show that the current in the ballistic 
regime is significantly larger than in the diffusive one. 

The spin current, in turn, is shown to have two com¬ 
ponents. One of them is related to the thermally induced 
spin polarization of electrons transferred from one reser¬ 


voir to the other, while the other one reveals the spin 
Nernst effect, i.e. the thermally induced spin Hall effect. 
We have also calculated the heat transferred by ballistic 
electrons from the reservoir of higher temperature to the 
reservoir of lower temperature. 
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